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The  Least  Mean  Square  algorithm  is  a method  of  extracting  a target 
signal  in  the  presence  of  interfering  noise  sources.  The  method  and  the 
computer  tests  for  simulated  conditions  utilizing  random  numbers,  are 
reviewed.  These  tests  indicate  that  the  method  is  feasible  for  laboratory 
conditions.  The  purpose  of  a sea  test  is  to  evaluate  the  method  for 
actual  ocean  conditions.  The  proposed  experiment  is  part  of  project 
PARKA  and  will  use  elements  of  the  Sea  Spider  array,  the  USNS  SANDS 
(AGOR-6),  and  the  Univac  1230  computer  system.  The  data  will  be  processed 
in  real  time.  The  experiment  and  computer  program  are  described. 


INTRODUCTION 
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In  this  memorandum,  the  proposed  experiment  of  adaptive  beamforming 
is  reviewed.  The  review  includes  a brief  description  of  adaptive 
beamforming,  the  sea  spider  array,  the  proposed  experiment  and  the 
computer  program. 


ADAPTIVE  BEAMFORMING 

We  will  follow  the  theory  developed  by  Widrow^)  and  Griffiths  ^ 
on  the  subject  of  adaptive  beamforming.  A tapped  delay  line  is  attached 
to  each  element  of  a K element  array.  There  are  L weights  in  each  tapped 
delay  line.  (See  Figure  1).  The  purpose  of  the  algorithm  is  to  adjust 
the  weights  such  that  there  is  a least  mean  square  (LMS)  difference, 
between  the  output  y and  the  expected  signal.  The  problem  was  first 
solved  by  Wiener  and  the  final  answer  was  in  terms  of  an  inverse  matrix. 
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Wldrov  and  Griffiths  obtained  an  approximate  iterative  solution.  The 
«uHr>  attractive  features  of  the  solution  ie  that  the  noise  field  need 
not  be  known  before  hand  and  that  it  is  not  neoessary  to  invert  a 
matrix.  In  our  simulation,  the  power  spectrum  of  the  signal  is  assumed 
to  be  known.  In  obtaining  the  IMS  difference,  the  noise  tends  to  be 
minimized  and  the  pattern  of  the  array  tends  to  have  notches  where  the 
noise  interferences  are  strongest. 

The  adaptive  beamforming  test  has  been  successfully  simulated  by 
the  use  of  random  number  data  by  students  at  Stanford  University,  The 
purpose  of  this  test  is  to  examine  adaptive  beamformer  performance  under 
actual  ocean  conditions. 
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SEA  SPIDER  ARRAY 

The  Sea  Spider  array  is  approximately  350  miles  north  of  Hawaii. 

It  is  a 3 legged  array,  each  leg  being  45°  from  the  vertical.  The 
array  is  buoyed  up  by  an  ellipsoidal  subsurface  float.  In  addition, 
over  3000  hollow  glass  spheres  are  mounted  on  the  three  legs,  making 
the  legs  neutrally  buoyant.  The  tension  in  the  3 cables  and  the 
neutrally  buoyant  elements  attached  to  them  combine  to  render  a stable 
and  rigid  support  for  the  array.  The  total  number  of  hydrophones  is  30 
and  there  are  10  elements  per  leg.  The  depth  of  the  water  at  the  array 
site  is  about  19,000  ft. 


PROPOSED  EXPERIMENT 


A linear  5 element  array  will  be  used  in  the  experiment.  The  four 
distances  between  nearest  neighboring  elements  are  21  ft,  52  ft, 

100  ft  and  171  ft. 

A 400  cps  projector  will  be  suspended  from  the  USNS  SANDS  (See 
Figure  2).  The  depth  of  the  projector  will  be  at  1000  ft.  The 
5 element  array  is  at  a depth  of  2500  ft.  The  ship  will  be  12,000  ft 
from  the  array  so  that  the  projector  is  in  the  Fraunhoffer  region  of 
the  array.  The  surface  reflected,  direct  and  bottom  reflected  rays  will 
arrive  at  28,7°,  37.9°  and  115.8°  respectively  from  the  broadside 
direction  of  the  array. 

The  Honeywell  projector  is  driven  by  a pseudo-random  noise  (PRN) 
signal  centered  at  500  Hz  in  a 200  Hz  band.  The  projector  transmitting 
response  curve  is  shown  in  Figure  3.  The  HX-90  has  a maximum  source 
level  of  100, 6db  (re  l^bar  at  1 yd).  Driving  the  projector  with  a PRN 
signal  reduces  the  source  level  to  91.1  db  and  at  a range  of  12,000  ft 
the  expected  acoustic  level  at  the  receiving  hydrophones  is  +19.1db. 
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A bio  ok  diagram  of  the  transmitting  system  is  shown  in  Figure  4. 

A PRN  signal  that  can  either  be  continuous  jar  gated  is  filtered  through 
a Buttervorth  filter  that  has  a 200  Ha  bandwidth  centered  at  500  Ha. 

The  signal  is  then  amplified  and  transmitted.  The  CML  amplifier  has  a 
maximum  output  power  of  5 KV  and  with  a 40%  projector  efficiency,  the 
electrical  power  is  more  than  adequate  to  drive  the  projector  at  full 
power  if  required. 

A block  diagram  of  the  receiving  equipment  is  shown  in  Figure  5. 

The  signal  reoeived  at  the  hydrophone  is  amplified  and  the  information 
is  telemetered  back  to  the  ship  where  it  is  filtered  through  identical 
Buttervorth  filters  and  reoorded  on  magnetic  tape  in  analog  form.  The 
analog  signals  also  go  to  a 12  bit  A/D  converter  and  these  ere  processed 
on  the  1230  Dhivac  computer.  The  digitised  signals  are  also  reoorded 
on  Univao  1540  magnetic  tape.  Baoh  execution  cf  the  IMS  Pattern 
Program  results  in  three  records  and  an  end  of  file  being  written  on  the 
Univao  1540  magnetic  tape  units.  All  records  are  reoorded  at  800  BP I, 
Odd  Parity,  Modulus  6,  and  in  Bioetal  Format.  The  order  of  the  three 
records  are  shown  in  Figure  6.  The  details  of  the  format  for  each  of 
these  records  is  shown  in  Figure  7 and  8. 


COMPUTER  PROGRAM 

The  updating  of  the  weights  is  obtained  from  equation  (l). 

W(i)  - Wti ) +t>Zzs{i)  - y X(if 7 
y - xu)  ^ 

C-  I 

where 

x(i)  are  the  KL  inputs 
y is  the  output 

pa(i)  is  the  signal  correlation 
U is  an  arbitrary  constant 
w(i)  are  the  KL  weights 

The.  inputs  x(i)  are  hard  clipped.  The  CS-1  program  in  fixed  point 
arithmetic  for  equation  (1)  is  given  in  appendix  A. 

The  beam  pattern  P(0)  for  a given  set  of  weights  is  given  by 
equation  (2).  See  reference  3 for  a full  explanation  of  this  equation. 
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The  subroutine  for  equation  2 la  given  in  Appendix  B. 


COMPUTER  RESULTS 

Random  Gaussian  numbers  were  used  to  simulate  a plane  wave 
striking  the  array  at  an  angle  of  22°  from  broadside.  The  power 
spectrum  corresponding  to  the  time  waveforms  was  flat  with  a center 
frequency  of  250  ops  and  a bandwidth  of  300  ops.  The  conventional 
pattern  of  the  equally  weighted  5 element  array  is  shown  in  Figure  9. 
Using  the  LMS  algorithm  and  a V =.01,  the  weights  of  the  tapped  delay 
line  were  updated  and  after  2000  iterations,  the  pattern  corresponding 
to  these  weights  was  computed.  The  result  is  shown  in  Figure  10. 

After  4000  iterations,  the  pattern  shown  in  Figure  11,  was  obtained. 
The  maximum  sensitivity  is  0 db  in  the  look  direction,  i.e.  the 
broadside  direction  in  this  oase.  A null  of  greater  than  -20db  is 


obtained  In  the  interference  direction 
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CONCLUSIONS 

A computer  program  using  CS-1  language  and  fixed  point  arithmetic 
has  been  written  and  is  operating  successfully  on  simulated  data. 

This  program  will  be  used  at  sea  to  test  the  applicability  of  adaptive 
beamforming  under  aotual  ocean  conditions. 
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Computer  Specialist 
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MAGNETIC  TAPE  FORMAT 
Figure  6 
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NOTES: 

1.  The  first  word  of  the  NOISE  record  is  "NOS"  left  justified  in  ASCII  oode. 

2.  The  first  word  of  the  SIGNAL  record  is  "SIG"  left  justified  in  ASCII  oode. 

3.  The  last  45  words  stored  in  the  upper  15  bits  of  X are  identical  to  the  first 
45  words  stored  in  the  lower  15  bits  of  X.  This  facilitated  programing. 

4.  X3,2008  is  interpreted  as  the  2009th  sample  on  channel  3. 

5.  Data  block  X is  comprised  of  5 channels  at  4009  samples  each  for  a total  of 
20,045  unique  data  words  in  A/0  code. 
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NTJSL  T«oh  Memo 
S«r  2211-162-69 


3£T *KHl  * 1 0* ( HHO ) l LHS i 
:»LH>SUl*,l*Tu*(HPO)  (SriNPI+abMl 
EwUi'Hl 

Y AH  Y *i\*F  t\OM*U*THHU«  IttOu 
SET*ld*Tw*K/3toOO 
SInCOS*Tl>*E*CwS 
SET*SU**|0*0 
StT*bO^2*To*0 
SET*SU»3«T0*0 

3E1 *COSC*To*COS/C 

WAkY*N*Fi<0M*0*THKU*3 

VANY*M*Fi<0M*M-1*THKW*4 

V AKY  * 1 *Ft<QM*U*THKU*9U 
SET  *A*TO# ( #P IM» I ) ) 

STK*A*ft(ftPMJ) 

SET*A*10*(*P(N.in 
STh*A#»*  (wPNl) 

SET  *CMNP*TO* ( »PHJ ) C WPNI ) ♦CMNP 
£N0*12 

POY**  < CNr^P ) *N  ( SMNP ) 

PUT*0*»  (CMnP) 

StT*T*TO*(QTlM»h) ) (COSO 
£NT*0*ft(oW) 

MULoh(l) 

KSH*AW*1'0XV  i»Y  * SC  IS  460 
LSH*au*UD»SC  IS  300 

HSH* A*30u * TRUNC M fc  WHOLE  OaMS  AND  EXTEND  SIGN 

KSH*Att«3*SC  IS  klJ 

STK*U*S(m1M) 

SlNCOS*AiM*S*SINl 

ENT*G*S(rO) 

MOL*fc(l) 

LSH*AO*l*U*SC  IS  300 

KSH*A*30u*THUNCATE  whole  BmMS  and  EXTEnO  sign 

HSH*ad*3»SC  IS  27 v 

STK*w**(a2M) 

3INC0S*A*M*E*C0Sl 
S£T*Al*TO*(PlbW)  U) 

SfcT#HPL*lO*(SiNl) (C0S1)/AI 

4F*T*fcG*U*rHtN*SfcT*KP0*T0*l 

3tT*SUtS2*T0*(hP0)  <SMwP>*30M2 

£Nu*m2 

END*N2 

SE1  *SUf-.2*T0*  UH  SON2 ) 

^AHY*P*FhOH***THHU*9U 
VAKY*M*FaOM*vKTHKU*4 
VAKt *N»F  nON*U*THKO*4 
iF*N*tu*H*THtU*6wT  U*0<' l 
VAKY*I*Fi\OKi*U*THKO*Gi>«P 
^ET*HPA*»0*P'fi 


APp«ndlx 


NUSL  Teoh  Memo  >■ 
Nov  2211-162-69 


b£1*A*TO*(WPiMfPPI)) 

St1*A*TO*UP(f«'I)) 

STft«A*NUPNl) 

SETwCMNH*T0*O»PMj)  (WPNI)+C*NP 
tNi)*U 

PUT*n<CMNP)*W<SMnP) 

PuT*U**(wMhP> 

SET*T*10»(PHL£L) 

SE1*T*T0*<DT<M»N))  ;COSC)«-T 
£NT*«*w(bW) 

MUL*W(H 

HSH*AQel  *01 V tiT  2 SC  IS  160 
LSH*AW*l20»SL  IS  300 

KSh*A*«30u'  TRUNCATE  WHOLE  SAMS  ANL  EXTEND  SJGH 

KSH*Au*J*SC  IS  27 0 

STK*V)*W(m1M) 

SINC0S*A1M*S*SIN1 

ENT*0*W<F0) 

HOL*W(T ) 

LSh*AL*UO»SC  IS  3U0 

HSH«A*^OU»TKONCATt  WHOLE  SAMS  AND  EXTEND  SIGN 

HSH*mW*3*SC  IS  27 u 

STR*0*w(A2M) 

SINCOS*AcM*E*COSj. 

SfcT*AI*Tw*(PIbW)CT>  * 

SET *KPu*T0* ( SlNi ) (vOSi)/Al 
IF*T*tu*0*THENwStT*RPOwTO*I 
SET*SUm3*T0*(KPD)  (SMNP) ♦aUM3 
ON  NO-OP 

tNO*NS 
ENUeHS 
EimU*P3 

SET *SUh3* TO*  ( 2 ) (SUM3) 

SE I *SUM*  I OeSUMl^SOHd+SUMJ 
PUT*W (SUi*»)  *W (POW+oI ) 

INCREMENT 

enlwio 

£NT*t>l*0 

LL*A 

PM  EnT*u*K  (POh+uI  ) 

COft*b*A*  i MOHt 
SIHeL*#» 

0Siv»U*  HuttO 
JPePw 

a Jh*A*n(nArt) 

Eh!  ♦li*cuL' 

PNA  «JP*cFaNuCH 

l»J!*UJ*P>’<1 
l.ut  Ah*c‘*w*HL«L» 
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NTJ3L  T«oh  Hmmo 
No.  2211-169-69 


-10 

NOKHA 


NltfW 


FOKM-TLXr*PLAb#l*btG  -SO  -4  s -4U  -35  -30  -25  -20 

5 0 OtClBtUS  wEl&rtT* 

LNT *A*PLAB 
KJPwftONRoE 
ENT*0*MPOw*ttl) 

Cl*  a 

LSH*AQ*lSO 
U1V**(NA  a) 

CONLOG IT* 'OUTPUT  IS  SCALED  UD 
STk*0*W (lPOm+01 ) 

STR*0*M1LP0M) »SC  U 1*L 
STh*ttl*MK) 

CLEAR *Z4u*PLAb 
FORM-Dt C *PL Ab*6 1 U* TLPO# 

ENf*A*49U 

SUB*A*BUAPOS 

UP*NW*f 

£NT*A*W (nT+B1) 

STR*A*MPRWT) 

F0f.M-DtC*PLAd*7lD*PRW  T 
FORM-DEC *PLAd* 1 *K 
ENT  *0*  W ( DOTS ) 

KPT*100*aDV 

STR*<i.*N<PLAB«H) 

EN7*Q*4l 

LSH*«*24l 

STR*&*W(PLAB*UO) 

ENI*0*M TLPON I • ALWAYS  NtoATlVE 

SUU*G*4000' ROUND 

RSH*0*l2u 

AQO*G*bOu*«POS 

ENT*U*b*SKIP 

AOU«u*b 

CL*A 

UIV*b 

STK*(j*lU64) 

Ei*T*05*A 
tNT*4**  (mASK) 

Cl*a 

KPT#u5 

K5h*Aw*0 

£KT*A**(r'LAb+t»4> 

SfcL*5U*MP0iUTX) 

STR***W(H.AB+t»4) 

ti4T*A*PLnU 

«JP*hOI<Ku£ 

USis*d1*1o00 

JP*NORMA 

RETURN 

LWL-pHuC*PATTEKL 
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Apponilx  B 


DEPARTMENT  OF  THE  NAVY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 

ARLINGTON,  VA  22217-5660  IN  REPLY  REFER  TO 

5510/1 
Ser  93/160 
10  Mar  99 


From:  Chief  of  Naval  Research 

To:  Commander,  Naval  Meteorology  and  Oceanography  Command 

1020  Balch  Boulevard 
Stennis  Space  Center  MS  39529-5005 

Subj : DECLASSIFICATION  OF  PARKA  I AND  PARKA  II  REPORTS 

Ref:  (a)  CNMOC  ltr  3140  Ser  5/1 10  of  12  Aug  97 

Enel:  (1)  Listing  of  Known  Classified  PARKA  Reports 

1 . In  response  to  reference  (a),  the  Chief  of  Naval  Operations  (N874)  has  reviewed  a number  of 
Pacific  Acoustic  Research  Kaneohe- Alaska  (PARKA)  Experiment  documents  and  has 
determined  that  all  PARKA  I and  PARKA  II  reports  may  be  declassified  and  marked  as 
follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  Chief  of  Naval  Research 
letter  Ser  93/160,  10  Mar  99. 

DISTRIBUTION  STATEMENT  A:  Approved  for  public  release.  Distribution  is 

unlimited. 

2.  Enclosure  (1)  is  a listing  of  known  classified  PARKA  reports.  The  marking  on  those 
documents  should  be  changed  as  noted  in  paragraph  1 above.  When  other  PARKA  I and 
PARKA  II  reports  are  identified,  their  markings  should  be  changed  and  a copy  of  the  title 
page  and  a notation  of  how  many  pages  the  document  contained  should  be  provided  to  Chief 
of  Naval  Research  (ONR  93),  800  N.  Quincy  Street,  Arlington,  VA  22217-5660.  This  will 
enable  me  to  maintain  a master  list  of  downgraded  PARKA  reports. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


PEGGY  LAMBERT 
By  direction 


Copy  to: 

NUWC  Newport  Technical  Library  (Code  5441) 
NRL  Washington  (Mary  Templeman,  Code  5227) 
NRL  SSC  (Roger  Swanton,  Code  7031) 
i/DTIC  (Bill  Bush,  DTIC-OCQ) 


LISTING  OF  KNOWN  CLASSIFIED  PARKA  REPORTS 


Operation  Plan,  Pacific  Acoustic  Research  Kaneohe-Alaska  PARKA  Experiment,  Undated,  ONR,  48 
pages 

(NUSC  NL  Accession  # 49531) 

Fleet  Research  Project  109  PARKA  II,  Undated,  COMAS WFORPAC-OPORD-303-69,  Antisubmarine 
Warfare  Force,  Pacific  Fleet,  Unknown  # of  pages 
(NUSC  NL  Accession  # 093561) 

Preliminary  Operation  Plan  Pacific  Acoustic  Research  Kaneohe-Alaska  PARKA  Experiment,  June 
1968,  ONR,  Unknown  # of  pages 
(NUSC  NL  Accession  # 023063) 

LRAPP  Briefing  Report  on  the  PARKA  Series,  May  1969,  MC  Report  001,  Maury  Center  for  Ocean 
Science  (ONR),  20  pages 
(NUSC  NL  Accession  # 023375) 

Bathythermograph  Traces  from  PARKA,  20  May  1969,  NUSL-TM-221 3-1 18-69,  7 pages 
(DTIC  # B952  259) 

Bathymetric  Strip  Charts  in  the  North  Pacific  Ocean  for  Project  PARKA  II,  20  June  1969,  Naval 
Oceanographic  Office,  Unknown  # of  pages 
(NUSC  NL  Accession  # 05 1 659) 

PARKA  II  Experiment  Utilizing  Sea  Spider  ONR  Scientific  Plan  2-69,  26  June  1969,  MC-PLAN-01, 
172  pages 

(DTIC  # B020  846) 

PARKA  I - Acoustic  Processing  and  Results,  28  July  1969,  USL  Technical  Memorandum  No.  2210- 
015-69,  NUSC  New  London,  115  pages 

(NUSC  NL  Accession  # 202993-001)  (NRL  SSC  Accession  # 85009134) 

A Scheduled  At-Sea  Simulation  of  Adaptive  Beamforming,  19  September  1969,  NUSL-TM-221 1-162- 
69,  23  pages 
(DTIC  # B026  991) f 

Biological  Data  Collected  on  the  PARKA  I Transit,  23  October  1969,  NUSL-TM-221 3-262-69,  15 
pages 

(DTIC  # B952  263) 

PARKA  I Experiment,  November  1969,  MC  Report  003,  Volume  1,  Maury  Center  for  Ocean  Science 
(ONR),  84  pages 

(NRL  Accession  # 466930)  (NRL  SSC  Accession  # 85004881)  (DTIC  # 506  209) 


ENCLOSURE  (1)  TO  CNR  LTR 
5510/1  SER  93/160  DTD 
10  MAR  99 


